Phase transition, spin-charge separation, and spin filtering in a quantum dot 



o 
o 

(N 
u 

< 



Michael Pustilnik^ and Laszlo Borda^ 

'School of Physics, Georgia Institute of Technology, Atlanta, GA 30332 
^Research Group "Theory of Condensed Matter" of the Hungarian Academy of Sciences, 
Institute of Physics, TU Budapest, Budafoki ut 8., H-1521, Hungary 

We consider low temperature transport through a lateral quantum dot asymmetrically coupled 
to two conducting leads, and tuned to the mixed-valence region separating two adjacent Coulomb 
blockade valleys with spin S = 1/2 and S = 1 on the dot. We demonstrate that this system exhibits 
a quantum phase transition driven by the gate voltage. In the vicinity of the transition the spin on 
the dot is quantized, even though the fluctuations of charge are strong. The spin-charge separation 
leads to an unusual Fano-like dependence of the conductance on the gate voltage and to an almost 
perfect spin polarization of the current through the dot in the presence of a magnetic field. 
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In a single-electron transistor setup the number of 
electrons in a quantum dot is controlled by the po- 
tential on the capacitively coupled gate electrode. At 
low temperature N is close to an integer at almost any 
gate voltage, except narrow mixed-valence regions, where 
adding a single electron to the dot does not lead to a 
large penalty in electrostatic energy. The distance be- 
tween these regions sets the scale for the dependence of 
measurable quantities on the gate voltage, which makes 
it convenient to use a dimensionless parameter A'o, the 
gate voltage normalized by this scale. In terms of A'o, the 
mixed-valence regions are narrow intervals of the width 



Ajv T/Ec < 1, 



(1) 



about half- integer values of iVo . Here Ec is the charg- 
ing energy and T is the tunneling-induced width of single- 
particle energy levels in the dot. 

In a typical experiment a dot is connected via tunneling 
jimctions to two massive electrodes jlj. At temperatures 
in the range F < T <C Ec, the conductance G is sup- 
pressed outside the mixed- valence regions, resulting in a 
quasiperiodic sequence of well-defined Coulomb blockade 
peaks in the dependence G{Nq) (3,01 • When T is further 
lowered, G{No) changes dramatically due to the onset 
of the Kondo effect At T ^ pairs of adjacent 

Coulomb blockade peaks merge to form broad maxima at 
N « odd integer, separated by smooth crossovers from 
the minima at iV ks even integer. Although G{No) at 
T > r is very different from that at T ^ F, in the mixed- 
valence regions both functions are featureless 0, 1^ ■ 

The evolution of G towards its low-temperature limit 
can be rather complicated. Indeed, GaAs quantum dots 
with odd N usually have spin S'=l/2 0,0. In this 
case the dependence G{T) is characterized by a single 
energy scale, the Kondo temperature Tr-; G(T) increases 
monotonically with the decrease of T at T ^ F 2] . How- 
ever, dots with even N often have spin 5 = 1 rather 
than zero 0, 0) 0- Kondo effect then occurs in two 
stages, controlled by two different energy scales, Tk and 
Tjf < Tr- 0. The resuhing G{T) is not monotonic: G 



first raises, and then drops again when T is lowered 
The dependence of the conductance on the Zeeman en- 
ergy B of an applied magnetic field is also non-monotonic 
and is characterized by the same scales Tk and T^ 0- 

The values of Tk and T^ depend on A^o and their 
ratio, in general, is not tunable. A notable exception 
occurs when the conductances of the dot-lead junctions 
are very different, i.e. when the width F = Fi+ Tji is 
dominated by the contribution from the lead with the 
stronger couphng to the dot, say, F^ » Fj^. (Note that 
conductances of the junctions are easily tunable in lateral 
quantum dot systems such as those studied in 0). It 
can be shown |^ that in this limit <C Tk for all Nq. 
In particular, in the vicinity of the mixed-valence region 
7k ~ Ffl, while Tk T l T'j^ ■ Accordingly, the second 
stage of the Kondo effect will not develop if 



Ffl<Cmax{r,B}<FL. 



(2) 



In this paper we show that under these conditions the 
conductance in the mixed-valence region between the 
Coulomb blockade valleys with 5=1/2 and 5 = 1 on the 
dot varies with A'o on the scale which is parametrically 
small compared with A^v, in striking difference with the 
conventional smooth dependence described above. 

The dependence of the conductance on i? at i? ^ T is 
qualitatively similar to its dependence on T at T ^ i? |3| ■ 
Since _B-dependence is much easier to understand, we 
concentrate here on the limit T — » (the effect of a finite 
T is briefly discussed towards the end of the paper). 

The first inequality in Eq. ^ allows one to take into 
account the coupling to the right lead in the lowest non- 
vanishing order in Fb/F/,. The conductance at any finite 
B is then given by 



3 = 5T+3i> gs^sin^Ss. 



(3) 



Here g — G/Gq is the conductance normalized by Go ~ 
(e^//i)Ffl/Fi, the largest value conductance per spin can 
reach for strongly asymmetric coupling to the leads; gs is 
the conductance (in units of Go) for electrons with spin 



2 



s, and 5s is scattering phase shift at the Fermi level for 
electrons with spin s in the left lead. 

In order to calculate the phase shifts, we model a quan- 
tum dot coupled to a single lead by the Hamiltonian 



H = Ho + Ht+Hd. 



(4) 



The first term here describes the electrons in the lead. 
For a lateral quantum dot it is sufficient to take into 
account only a single propagating mode 0, Q , 



ks 



(5) 



and the spectra can be linearized near the Fermi level, 
which corresponds to a constant density of states v. 

The second term in Eq. @ describes the tunneling 
coupling between the dot and the lead. 



H.c. 



(6) 



nks 



In the following we set i„ = i, so that all levels in the dot 
have the same width F ~ irvt^. This assumption is not 
essential for the validity of the following consideration. 
The last term in Eq. describes an isolated dot. 



„,d„, + Ec{N- NoY - EsS' - BS,. (7) 



Here N = J^ns^isdns and S = ^J2nss' 4is^ss'd^^, are 
operators of the total number of electrons on the dot, and 
of the dot's spin, respectively ((T = (cti:, CTy, tiz) are Paufi 
matrices). For a typical dot the parameters 6E (mean 
single-particle level spacing) , £'5 (exchange energy) , and 
Ec (charging energy) satisfy Es <^ 5E <^ Eq [l] ■ 

An isolated dot with even N will have 5* = 1 in the 
ground state if the spacing e between the two single- 
particle levels closest to the Fermi level is anomalously 
small, 2Es — e > 0, IE ■ 1^°^ simplicity, we assume 
here that 2Es — e F. Although this simplification 
imposes a stronger restriction on F then F ^ 6E (this 
inequality justifies the tunneling Hamiltonian description 
of the dot-lead junction 2J), it does not affect the results. 

For the model Q-Q the phase shifts are given by 

6^ = {^/2){N + M), 5 1 = [■k/2){N - Af), (8) 

where N — (N) is the number of electrons in the dot, 
and M = 2{Sz) is the dot's magnetization 0. We start 
with A'o outside the mixed-valence region, 

Aat < |iVo -^ol -Cl, A^o - 1/2 = odd integer, (9) 

Here N No±l/2 is close to an integer. The tunncling- 
induced virtual transitions to states with "wrong" N can 
be "integrated out" with the help of the Schrieffer- Wolff 



transformation, yielding an effective Kondo Hamiltonian 



H = Ho + Vp + J{s-S) - BS, 



(10) 



where p = J^kk's -^L^k's and s = i T,kk'ss' ^L^^s'^k's' 
are operators describing the local particle and spin densi- 
ties of conduction electrons. The operator S in Eq. HlOfl 
is a projection of S [see Eq. on the ground state 
multiplct of an isolated dot with fixed integer TV. The 
reduction of the microscopic model Q-lO to the Kondo 
Hamiltonian Hl()|l is valid only when A^o is outside the 
mixed-valence region and at sufficiently low energies, 



161 <^ = min{djv, 2Ec\No-No\], 



(11) 



where (In = SE {(In ^ 2Ec - e) for odd (even) N. The 
parameters V and J in Eq. can be estimated as 0,0 



fVr^ An{No~ No) \ iyJr^\iyV\. 



(12) 



(It should be noted that V and J are subject to strong 
mesoscopic fluctuations; the order-of-magnitude estimate 
l|12|l is sufficient for our purpose). 

The potential scattering term in Eq. is responsible 
for the deviations SN of the dot's occupation from the 
corresponding integer values Nq ±1/2, 

SN = N~{Na±l/2) w -2iyV - AN{Na-Noy\ (13) 

Note that |i5iV| is finite and increases with approach to 
the mixed-valence region. Also note that a weak mag- 
netic field B < F does not affect N{No). 

On the contrary, M{Nq) depends strongly on B. In- 
deed, M{B) for a given A^o is controlled by the Kondo 
temperature Tk{No), which can be estimated from 



ln{D/TK) {vjy 



(14) 



Accordingly, Tk F at |7Vo A/^oj Aat, and decreases 
exponentially with the increase of the distance to A^o- A 
fixed field B is large compared to Tk at | A'^o — -^o | ^ Ab, 
where As is the distance between A'o and A'^o at which 
B ^Tk{Nq). In this regime 



M/Afo = 1 - 2\n{B/TK) 



(15) 



where Mq = 1(2) for A^o < A^o (A^o > A^o)- In the oppo- 
site limit As > jA'^o - Nq\ > Aat (note that As > An 
for _B -C F) the system is in the strong coupling regime 
B ^ Tk- Here M{Nq) depends strongly on the parity 
of A^. Indeed, S in Eq. H1U|I is spin-1/2 operator for odd 
A^ (i.e. for Nq — Nq ^ An), and spin-1 operator for 
even A^ {Nq — Nq ^ An)- This difference is crucial. 
An antiferromagnetic local exchange interaction with a 
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single species of itinerant electrons suffices to completely 
screen S = 1/2 magnetic impurity, ther eby forming a 
singlet (non-degenerate) ground state [ifl ll^ • In 

this case the approach to the low-energy fixed point is 
Fermi- liquid-like |llj . and 



M - B/Tk- 



(16) 



On the contrary, for S" = 1 only half of the impurity's 
in is screened, and the ground state is a doublet |lll 
The low-energy physics is then described by the 
ferromagnetic exchange of the conduction electrons with 

the remaining spin 5' = 1/2 |ilE3. 

and 



M = 1 



21n(rK/B) 



(17) 



Although the above results were obtained for A^o out- 
side the mixed- valence region |iVo — iVo| < Ajv, some 
conclusions regarding this region can be drown as well. 
Indeed, since 6N in Eq. (|13|l is finite, it is plausible that 
N varies continuously with Nq, as sketched in Fig.^a). 
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FIG. 1: (a) Number of electrons in the dot A'' differs appre- 
ciably from an integer in a narrow mixed-valence region of the 
width Ajv. (b) At _B ^ F, the width Am of the crossover 
region in the dependence of the magnetization M on the gate 
voltage No is small compared to Ajv. (c) Spin-resolved con- 
ductances in the crossover region at T ^ B ^ F. (d) The 
total conductance at max{_B, T} ^ F. 

The dependence M(No) is more complicated. Con- 
sider the limit B +0. In this limit AI is determined 
solely by the ground state degeneracy. Since the degen- 
eracy can not change continuously, the system must go 
through a quantum phase transition (QPT) at a certain 



value of A^n 



Ao*. 



As shown above, the ground state 



is either a singlet or a doublet when the charge fluctu- 
ations are weak. Therefore, the transition must occur 
within the mixed- valence region, i.e. |A^q — Aq | < A^r. 
The QPT manifests itself in a singular dependence of the 
magnetization M on the gate voltage. 



lim M = 6'(An 

3^+0 



A*). 



(18) 



Note that in the vicinity of the transition the spin is quan- 
tized, even though the fluctuations of charge are very 



strong, A'*= N{Nq)k half-integer (unlike in the case of 
transitions that occur at a fixed integer A^ 0,E|). 

Any finite field lifts the degeneracy of the ground 
state. QPT then turns to a crossover, and the sharp 
step in the dependence of M{No) is smeared. The 
crossover takes place in a narrow interval of gate volt- 
ages |Aq — A"q I < Am- We expect that at a sufficiently 
low field the crossover width Am remains to be small 
compared to An, see Fig. ^b). In order to estimate 
Am, we now construct an effective Hamiltonian -ffqpT 
for the vicinity of the transition. 

Such Hamiltonian should be applicable at low energies 
{B,T < F) and for A^o in the range \Nq - N*\ < An, 
which includes the crossover region. At these energies 
and gate voltages the number of electrons in the dot is 
approximately constant, A « A*, while M(Ao) changes 
rapidly. It is therefore plausible that Hqpt acts only 
on the spin degrees of freedom (spin-charge separation) . 
At energies below F half of the dot's spin when it is in 
the triplet state is already screened. The simplest pos- 
sible model accounting for the interaction of the (still 
unscreened) spin-1/2 with electrons in the narrow strip 
of energies \£,k\ ^ P reads 



H, 



QPT 



Hi 



BS"" 



(19) 



Here Hq and s' [cf. Hq and s in Eq. H1(J|) ] are defined in 
terms of the operators i/'L acting in the basis of single- 
particle states that incorporate an extra scattering phase 
shift S* = ttN*/2. 

For Hqpt to describe the change of the ground state 
symmetry at A'g = A^g, the exc hang e constant J' must 
change its sign at this point [10, lU ll^ ■ Assuming the 
dependence J'(A^o) to be analytical, we can write 



,yJ'{No)^A^\N*-N,). 



(20) 



The coefficient in Eq. (|20|l has been chosen in such a 
way that vJ' ^ 1 for Aq — A'g ^ An- This ensures the 
continuity of M(Ao) throughout the singlet side of the 
transition Nq < Nq. 

A comment on the status of Eqs. (|19|l . (|20|l is in or- 
der here. The effective low-energy Hamiltonian -ffqpT is 
in the same relation to the original microscopic model 
Q-Q) as, e.g., the effective Fermi-liquid description of 
strong coupling regime is to the Kondo model. As in 
the latter case, the applicability of iJqpT can be verified 
by comparing the predictions of the two models. 

The magnetization for the model Ijigil . (|20|l is obtained 
using the standard scaling arguments 10]. Very close to 
the transition (when \i^J'\ln{T/B) ^ 1) the first order 
perturbation theory in vJ' <^ 1 yields 



M - 1 « -I. j'/2 ^ {Nq - a*) . 



(21) 



On the doublet side M(Ao) slowly increases with the 
distance to the transition, saturating at 



M =1+ [21n(F/S)]" 



(22) 



Note that Eq. matches Eq. ((T7jl at the border of the 
mixed-valence region. 

On the singlet side of the transition {Nq < Nq) the 
magnetization is given by Eq. H15|) (with Mg — 1) for 
B :^Tk and by Eq. ITHJl for B ^Tk, where the Kondo 
temperature Tr-(A^o) satisfies 



ln(r/T;^) = (lyfy 



An{N*-N, 



0) 



(23) 



Tk increases with the distance to the QPT from Tk — 
at A^o = -^0 to Tk ~ F at the border of the mixed valence 
region, where it matches Eq. (|14|l . 

As iVo is tuned through the mixed-valence region M 
grows monotonically from M ^ B /T <^ 1 to the value 
given by Eq. (|22() . The increase takes place mainly in a 
narrow interval on the singlet side of the transition where 
B < Tk{No). Eq. ^ then yields the estimate 



HT/B) 



(24) 



The evolution of the phase shifts with A^o can now be 
deduced from Eq. |(SJ. To the left of the crossover [see 
Fig. n^b)] both phase shifts arc given by Sg ~ ttN* /2 
with N* Ki No, see Eq. ©. As iVo is tuned through 
the crossover, raises, while Si drops by approximately 
7r/2. Therefore the phase shifts necessarily pass through, 
respectively, the anti-resonance i5| = 0(mod7r) and res- 
onance Si — 7r/2 (mod tt). Hence, within the crossover 
region the conductances (O satisfy g^/gi ^ 1, and 
vanishes identically at some value of Nq- In other words, 
the system acts as a perfect spin filter. 

Details of the dependencies gs{No) are sensitive to 
the dot's occupation at the transition N* . While N* 
is close to a half-integer, it's precise value is obviously 
non-universal. For example, iV* depends on the values 
of tn for all n in Eq. ©. In Fig. QJc) we sketch gs(A^o) 
for < a ^ 1, where a = N* — Nq. The dependence of 
the total conductance g on A^'o in this case has a charac- 
teristic Fano-like shape, see Fig. ^d). 

In order to verify the applicability of the effective 
Hamiltonian H19|l . we performed extensive numerical 
renormalization group (NRG) |il3i] simulations. For this 
purpose, we truncated the dot's Hamiltonian to that 
of a two-level system with 



ne/2, n — ±1. 



(25) 



The NRG data, see Fig. |21 are indeed in an excellent 
agreement with the behavior expected from Eq. (|19|l . 
The sharpening of the step in the dependence M{No) 
with the decrease of B, obvious in Fig.|2a), is described 
very well by Am/ Aat = a[ln"^(r/B)-h foln"^(r/B)] with 
a = 3.0 and b — 9.5; at low field this agrees with Eq. (Pljl . 
Here we defined Am as the distance in Nq between the 
points where M = 0.5 and 1, and A^f as the distance be- 
tween the points in Fig.|2Ib) where N = 1.25 and 1.75. 
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FIG. 2: Results of NRG simulations of the model l|H-lf7jl. l|^ 
with r/Ec = 0.05, Es/Ec = 0.16, and e/Ec = 0.1. 
(a) Magnetization M{No) at different B. At the transition 
M — 1 independently of _B, in agreement with Eq. I|21|l . (b) 
Dot's occupation N(No). (c) Conductance g{No) at B/T = 
2.5 X 10"". Note that in this case a = N'-3/2 « -0.25 < 0, 
hence the difference with Fig.^d). 



So far, we considered the conductance at T = 0. The 
above results are valid as long as T ^ B; corrections 
to gs in this case are of the order of {T/B}'^, and the 
spin- filtering property remains intact: mm{g^ / gi} ~ 
(r/B)2 < 1. At T > S the field has a negligible ef- 
fect. The dependence g{No) in this limit is very similar 
to that at B > T, see Fig. IHd) andE^c), with T re- 
placing B in the crossover width Eq. H24|l . This peculiar 
dependence will be observable already at moderately low 
temperatures T < T (note that the observability of the 
conventional Kondo effect requires T < minjrR-} <C F). 

To conclude, we studied a lateral quantum dot asym- 
metrically coupled to two conducting leads, and tuned to 
the mixed-valence region between the Coulomb blockade 
valleys with = 1/2 and S" = 1 on the dot. This regime 
can be realized in devices such as those studied in ^. 
We predict that, contrary to naive expectations, the con- 
ductance varies with the gate voltage on the scale which 
is parametrically small compared with the width of the 
mixed-valence region. 
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